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Electron transfer coupled to a collective vibronic degree of freedom is studied in strongly condensed 
phase and at lower temperatures where quantum fluctuations are essential. Based on an exact 
representation of the reduced density matrix of the electronic+reaction coordinate compound in 
terms of path integrals, recent findings on the overdamped limit in quantum dissipative systems 
are employed. This allows to give for the first time a consistent generalization of the well-known 
Zusman equations to the quantum domain. Detailed conditions for the range of validity are specified. 
Using the Wigner transform these results are also extended to the quantum dynamics in full phase 
space. As an important application electronic transfer rates are derived that comprise adiabatic 
and nonadiabatic processes in the low temperature regime including nuclear tunneling. Accurate 
agreement with precise quantum Monte Carlo data is observed. 



I. INTRODUCTION 

Electron transfer (ET) in complex environments is ubiquitous in physics, chemistry, and biology, the most prominent 
example being photosynthesis 0, • In the last decade with the advent of pump-probe techniques the detection of 
fast electron transfer dynamics in molecular systems became accessible and meanwhile even coherent vibronic motion 
interacting with transfer processes has been observed in a variety of condensed phase structures |^|-|14|. Experi- 
mentally, the transfer dynamics depends very sensitively on a bunch of parameters such as temperature, electronic 
couplings, structure of the environmental degrees of freedom, to name but a few. 

Despite this complexity theoretical foundations to describe these systems have already been laid in the 80ties with 
the pioneering work by Marcus et al. [lj . It turned out that qualitatively two dynamical domains must be distinguished: 
For vibronic dynamics fast compared to the bare ET, characteristic for long-distance ET, one speaks of nonadiabatic 
reactions, in the opposite case of very slow vibronic motion, realized e.g. in mixed-valence compounds, the reaction 
is said to be adiabatic. Tremendous progress towards a deeper understanding of this picture has been made for 
the archetypical case of donor-acceptor systems, especially based on descriptions that include the coupling to the 
residual molecular degrees of freedom also at lower temperatures. The common strategy starts from system+reservoir 
models and considers the reduced dynamics of the system (relevant degrees of freedom) only In this context 

a standard treatment has been Redfield theory H^-jl^, but also path integral approaches have been exploited 
|l5|. the latter one particularly suitable for strongly condensed phase systems. Within these formulations transfer 
rates have been successfully calculated in the nonadiabatic regime over the whole temperature range by invoking 
golden rule techniques. Likewise, rate calculations for adiabatic reactions at sufficiently high temperatures exploited 
well-established tools for activated rare events as e.g. Kramers' flux over population method. These quite different 
approaches reflect the different physical processes that control the transfer, namely, in the former case the electronic 
coupling between diabatic donor/acceptor states, in the latter one the sluggish activated bath motion on the lower 
adiabatic surface. To find a formulation that interpolates, at least partially, between these distinct dynamical ranges 
has been one of the fundamental problems in ET theory. 

A crucial step forward into this direction has been first achieved by Zusman [2(| and then later within the path 
integral approach by Garg et al. [2l|. For sufficiently high temperatures and moderately fast to slow bath modes they 
derived equations of motion for the electronic dynamics interacting with a damped collective degree of freedom, coined 
the reaction coordinate (RC). In essence, a classical overdamped RC motion (Smoluchowski limit) is coupled to the 
Heisenberg equations of the electronic degree of freedom. These so-called Zusman equations (ZE) not only opened the 
door to study the relaxation dynamics of the compound in detail H^-H^, but also to derive an explicit expression for 
the relaxation rate [Hl20||. The Marcus/Zusman rate comprises to some extent nonadiabatic and adiabatic dynamics 
as well, and reduces to the high temperature golden rule result for small electronic coupling and to the corresponding 
activated Kramers rate for large coupling. A thorough discussion of the various rate expressions and their performance 
in comparison with precise quantum Monte Carlo results has been given in |26j |. 
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The only range where so far no analytical description has been available includes moderately fast to slow bath modes 
as in the ZE, but at low temperatures. For a RC this domain corresponds to overdamped quantum dynamics which 
has been analyzed only recently in pTl l28l| : There it has been shown that the position probability distribution of a RC 
obeys a Smoluchowski type of equation also at lower temperatures with a substantial impact of quantum fluctuations; 
for high temperatures the quantum Smoluchowski equation reduces to the known classical one. The basic idea of the 
present work is now to investigate to what extent these findings can be adapted to the ET problem and to arrive at a 
generalization of the ZE to the quantum domain. For that purpose the path integral representation must be utilized 
meaning here to start with the exact density matrix of the reduced dynamics for the electron+RC. Fortunately, the 
corresponding rather involved calculation has already been carried out in [23. In the sequel, we directly built on 
these achievements and derive in a kind of semiclassical analysis from the exact expression approximate generalized 
Zusman equations (GZE). One particular advantage of the approach is that the approximations are well controlled 
and detailed conditions for the range of validity of the GZE can be given. 

The article is organized as follows. We introduce the model and collect the main results for the reduced dynamics 
in Sec.|n| In Sec. IHII we start with a discussion of the quantum Smoluchowski limit to proceed with the derivation of 
the generalized Zusman equations from the exact path integral expression. An extension to the full RC phase space 
in terms of Wigner functions is achieved in Sec. II VI As an application we derive in Sec. an analytic expression for 
the ET rate which generalizes the classical Zusman result to lower temperatures and compares well with numerically 
exact quantum Monte Carlo data. Conclusions are given in Sec. IVII 



II. MODEL HAMILTONIAN AND REDUCED DENSITY MATRIX 



We consider intermolecular ET between localized electronic states, called donor and acceptor henceforth, in presence 
of a vibronic environment. Basically two models have been developed to describe the dynamics of this compound, 
namely, the spin-boson model [23] and the reaction coordinate (RC) model [23. Within the reduced subspace of the 
electronic dynamics both models are completely equivalent, but they differ in the representation of the environmental 
degrees of freedom and thus allow to focus on different aspects of the transfer process. In the spin-boson model the 
electronic two state system is linearly coupled to a bath of harmonic oscillators, while in the reaction coordinate picture 
the electronic system interacts with a collective bath degree of freedom (reaction coordinate) which is embedded in 
a harmonic environment. The latter modeling allows to study explicitly the combined dynamics of electronic and 
collective vibronic motion and particularly the influence of vibronic wave packet motion onto the transfer process. 



A. Electron Transfer Model 



The corresponding reaction coordinate Hamiltonian consists of four parts, namely, 

H = H EL + H RC + H! + H b (1) 

where 

he hA 

Hel = -—<Jz —o x (2) 

represents the bare electronic two-state system as an artificial spin-i system with donor |— ) and acceptor |+), a bias 
he, and an ET coupling hA which is assumed to be independent of the nuclear degrees of freedom. The second part 

Hrc = ^- + ^tV " C q<J z 3 

describes a harmonic collective vibronic degree of freedom (reaction coordinate, RC) linearly coupled to the electronic 
system and interacting with a heat bath of residual molecular degrees of freedom 

c 2 
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This way the environmental oscillators couple to the transfer system only indirectly via the reaction coordinate q 
(secondary bath). Of course, one could transform both q and {yj} to normal modes which leads to the spin-boson 
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model |21j,|31|. Here, however, we are interested in the combined dynamics of spin+reaction coordinate and thus need 
to calculate the time-dependent reduced densit y m atrix p acr i (q, q', t) obtained by integrating out the secondary bath. 
This can be done exactly as we have shown in [23] and we will summarize the main findings in the next section. 

Before we do so, let us address the initial condition for the time evolution which turns out to be a crucial point, 
especially in the strongly condensed phase limit. We imagine a donor state as being excited from some ground state 
(dark state) at time t = (see fig. 1). The excitation pulse is assumed to be sufficiently short so that there are 
no correlations between the electron and its surrounding at t = meaning that the total initial density matrix Wq 
factorizes into a direct product of a projector onto a pure donor state and a density operator po of RC and secondary 
bath, i.e., 



W 



-)<-| <8 Pa- 



(6) 



The generalization to other electronic initial states is straightforward (see below). For po two qualitatively different 
cases have to be distinguished. In the first one RC and secondary bath factorize corresponding to an uncorrelated 
initial state. This "preparation" is well-known from Feynman- Vernon theory [32| and is clearly reasonable if RC and 
bath are only weakly coupled. In case of strong coupling, however, it is more realistic to assume a correlated initial 
state 33] 



pl or {q-,q \{vn},{y' n }) = pp(q,{yn},q',{y' n }) 

with the equilibrium density matrix of RC and secondary bath in the donor state 



Zb 



RC 



(q,{yn}\exp 



-a h 



H, 



P 

2m 



co q 



\q',Wn}) 



(7) 



(8) 



at inverse temperature f3 — and with .Zb,rc the partition function of the bath+RC compound. Deviations 

from equilibrium in the RC are created by a preparation function A(g, q'). In fact, only starting with correlated initial 
states guarantees to regain from the exact quantum dynamics of the RC in the high temperature/strong damping 
limit the so-called Smoluchowski equation for its position probability distribution [34 . l35| . 



B. Reduced Dynamics 

Exploiting Feynman's path integral formalism for the time propagation of the reduced density 

p{t) = Tr B {exp{-iHt/h)W enp(iHt/h)} (9) 

the harmonic degrees of freedom of the secondary bath can be eliminated exactly j3|| |3tJ • This leads to the exact 
expression 

p tTft , / (q f ,q' f ,t) = r (t) ~ a ' Va f ^ 2V e *C^M-^M) ^ q > t . H [(T ']) . (i ) 

Jcr(0) = -1 Jcr'(0) = -1 

Here, forward and backward electronic path a(s),a'(s) connect in time t the initial points <j(0) = — l,o"'(0) = — f 
with <Jf,a'j: where each path is weighted with the bare electronic action 5el and the density matrix of the damped 
RC. The latter takes the form 

1 f°° 

/5(?/)9/,*;HH) = J dq t dq[ \(q t ,q' t ) J(g/, g}, i, ft, g-; [ct], [cr']) (II) 

where the propagating function J(-) is a threefold path integral (two in real time, one in imaginary time) over the 
RC degree of freedom only. The forward and backward real time paths q(s),q'(s) run in time t for given electronic 
paths a{s) 1 a'{s) from qi,q[ to fixed end points qf,q'f, while the imaginary time paths q(r) connect qi with q[ in the 
interval h(3. The contribution of each path is weighted according to its effective action Erc, i-e., 

^Rc[q,q, q';cr, a'] = iS RC [q] + Src [q; <j, a'] - Snc[q' ;&,<?'] + <f>[q,q,q'] ■ (12) 

It consists of the bare actions of the RC system according to Hrc © in real and imaginary time and the so-called 
influence functional describing the effective influence of the secondary bath onto the RC. The real time paths encode 
the dynamics of the RC system and the imaginary time paths specify its initial preparation, namely, 

p(qf,q' f , 0; <n, °[) = -^Hqf,q' f ) pp(qf,q'f) (13) 
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with pp{q : q') — (<z|Tre exp[— j3(H^c + Hb + Hi)]\q')\ a ^ = _i the reduced equilibrium of the RC in the donor. The 
normalization in 1|12|) and (|13fl is chosen such that for vanishing RC-bath coupling the influence functional reduces to 
1, i.e. Z = Zb,hc/Zb where Z-q is the part ition function of the secondary bath alone. We omit here the explicit form 
of the effective action and refer to p9l l33j for further details. Note that since electronic and RC-system are initially 
uncorrelated the Euclidian action <Src[?] in l|12l) is independent of the electronic paths. The influence of the harmonic 
bath onto the RC described by the influence functional <j)[q, q, q'\ is completely covered by the autocorrelation function 
of the fluctuating force ^CjUj [see (0] 







3C 



dio ^ cosh[a;(ft/?/2-*?)] 

- I{UJ) smh(cohf3/2) (M) 

where 9 = s — it, < s < t, < t < Kj3. All relevant bath properties are contained in the spectral density I(u>) of the 
bath oscillators which, as a classical quantity, can be obtained from molecular dynamics simulations. For real times 
K(s) = K'(s) + iK"(s) is directly related to the macroscopic damping kernel 

7(s) = — / cos(ws) (15) 

mJ Q it to 

via K"(s) = (m/2)d r y(s)/ds and K'(s) — * mj(s)/hp in the classical limit cj c ?L/3 <C 1 with a typical bath frequency u> c 
usually taken as the bath cut-off frequency. The damping kernel produces interactions along the RC paths that are 
non-local in time even in cases where u> c is very large and thus, classically, the corresponding bath memory time very 
short. This is due to quantum fluctuations that appear on a time scale Kf3 and are responsible for non-Markovian 
behavior at lower temperatures. 



C. Exact Reduced Density Matrix for the Reaction Coordinate 



For given electronic paths the dissipative dynamics of the RC looks like the dynamics of a damped harmonic 
oscillator in presence of force fields |33j. The crucial difference to the case of an external force, however, is the fact 
that electronic forward and backward paths are not identical, but rather are dynamical degrees of freedom as well. 
The corresponding somewhat involved calculation to obtain the reduced density of the RC has been presented in [2^ . 
Eventually, one ends up with an exact expression for the reduced dynamics of the RC 



P(xf,r f ,t; [x\,[v]) = 



1 



N(t) 



dxi dri \(xi,ri) exp 



-T, ma .(xf,rf,t,Xi,ri; [x], [17]) 



where for convenience we introduced sum and difference paths for the spin 

X =(a-a')/2, V =(a + a')/2 
and sum and difference coordinates for the RC 



(16) 



(17) 



x a =q a ~q' al r a = (q a + q' a )/2, a = i,f. (18) 

The exponent S ma is the effective action evaluated at the minimal action paths and N(t) is the normalization. The 
explicit form of both is not very illuminating (for details see |29j): E ma is a quadratic function of the RC-coordinates 
with time dependent coefficients that are functionals of %(s) and 77(a), while N(t) depends only on time. For given 
spin paths their time dependence can be completely expressed in terms of the real and imaginary parts of the 
autocorrelation function of a damped harmonic oscillator (q(t)q). 

The exact dynamics of the electronic+RC system in presence of the residual vibronic degrees of freedom is now 
given by (flT7f> together with l(TH|l . If one would integrate over the diagonal part p(0, rj,t; [x], [rj]) too by putting A = 1, 
one would recover the exact reduced path integral expression for the reduced spin-boson dynamics. Here, however, 
we are able to study the relaxation dynamics of the RC starting with nonequilibrium initial states explicitly. For 
that purpose what remains to be done, is an evaluation of the spin path integrals which is not possible analytically. 
The fundamental complication lies in the self-interactions, non-local in time, and mediated by the secondary bath. 
As a consequence, the dynamics of the density matrix p af!(T i(xf,rf,t) can in general not be determined from simple 
equations of motions so that one has to work within the path integral representation. Progress can be made by 
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applying either numerical methods like Quantum Monte Carlo techniques |26| or perturbative approaches like e.g. 
NIB A or golden rule [l^. As was shown by Zusman [2(j and then in detail by Garg et al. [2l| particularly the high 
temperature/strong damping limit allows for an effectively Markovian description. There, only the diagonal part 

P<r f ,*' f (qf,t) = Pa s ,a' s {Xf =0,r/ =qf,t) (19) 

is relevant and its time evolution is determined by a set of equations of motion where the bare electronic dynamics 
is coupled to a classical overdamped RC motion (Smoluchowski limit). Our first goal is to derive generalized Zusman 
equations which are also valid in the low temperature range. 



III. GENERALIZED ZUSMAN EQUATIONS 



While the strong friction limit is well-known in classical mechanics |38| , corresponding quantum systems have been 
analyzed only recently |27|. Hence, in this section we first summarize some main findings which then serve as a basis 
for the treatment of the ET system with overdamped RC. 



A. Classical and Quantum Smoluchowski Limit 



The strong friction limit in classical physics, known as the Smoluchowski limit, is characterized by a separation of 
time scales between equilibration of momentum and equilibration of position. This allows to adiabatically eliminate 
the momentum from the phase space Fokker-Planck equation and to gain a simple time evolution equation for the 
position distribution, the Smoluchowski equation |39j . In case of a harmonic oscillator with frequency ujq the basic 
condition then is 7 ^> ujq where 7 denotes a typical damping strength. This is defined as 

7 = 7(0) = fimM (20) 

with 7(w) the Laplace transform of 7(f) (jl5|l . For instance, in case of ohmic friction I(uj) — jmu> and also for 
the more realistic Drude damping (also called Debye spectral density) I(u>) — mjLuuj c / \uj 2 + u)%) one has 7 = 7. 
Since one is interested in the dynamics on long time scales where the position distribution equilibrates, only the low 
frequency modes of the residual degrees of freedom are relevant. Now, for overdamped quantum systems given a 
typical frequency luq of the bare system, e.g. its ground state frequency, by strong damping we mean 

T 11 

4, » ft/3,-,-. (21) 

^0 7 

In other words, the time scale separation known from the classical Smoluchowski limit is extended in the quantum 
domain to incorporate also the time scale for quantum fluctuations Kf3. Correspondingly, one considers the dynamics 
on a coarse grained time scale s 3> h(3, — in real time and r 3> — in imaginary time. The consequences are 
substantial: (i) The strong friction limit suppresses non-diagonal elements of the reduced density matrix during the 
time evolution. This simply reflects the fact that a quantum system tends to behave more classically, the stronger 
coherences are destroyed in presence of a heat bath, (ii) The real-time part of the damping kernel K (s) becomes local 
in time so that a time evolution equation of the form p(t) = C pit) with a time independent operator C may exist. For 
the quasi-ohmic case uj c 3> 7 considered in the sequel the corresponding range in parameter space covered by l|21|) is 
shown in fig. 2. It is well separated from the weak friction limit region and comprises temperatures from the classical 
"fh/3 <C 1 to the deep quantum domain 7 ft/3 S> 1. 

This way it has been shown in |2?| that for continuous systems moving in sufficiently smooth potentials V(q) the 
diagonal part of the reduced density matrix P(q,t) = p(q,q,t) obeys the equation of motion P(q,t) = £qs P(q,t) 
with 



C ^-^8q{ VM+ 0-q 
where ■ abbreviates d/dt, 1 stands for d/dq, and V c g — V + AV" /2. Here 



J,+W"(q) 

P 



(22) 



2 °° I 

A = hp § Vi + v^{v n ) (23) 



=1 
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with Matsubara frequencies v n — 2irn/hl3 measures typical quantum fluctuations in position space. In the high 
temperature range yhfl <C 1 these are negligible and one recovers the classical Smoluchowski operator Cqls — £qs (A = 
0). For low temperatures, however, quantum effects appear and for jhfi 3> 1 we find A « (7t/m7)rn(77i/3 /2tt) . In the 
specific case of a harmonic oscillator the diffusion coefficient in Cqs becomes independent of position and differs from 
the classical one just by a renormalized temperature 1//3 — > 1//3 + Ajtiwq. Moreover, one has (q 2 ) ~ (9 2 )ci + A. 

Indeed, the strong coupling treatment even extends to classical phase space so that the Wigner transform W(p, q, t) 
of the density matrix p(q,q',t) satisfies a quantum Fokker-Planck equation psf . As the quantum fluctuations in 
position are squeezed by friction, the quantum fluctuations in momentum are enhanced. Namely, for jhft 3> 1 they 
grow like (p 2 ) « (mhj / ir)ln(uj c /~f) . 

Based on these results we now want to derive from the exact density matrix p r7fr7 ^(xf,rf,t) in (|I0|l approximate 
equations of motion for its diagonal elements in case of an overdamped RC. We know that in the classical realm 
7?l/3 <C 1 this leads to the so-called Zusman equations (ZE). A naive expectation would be just to replace the classical 
operator for the RC motion C c \ by its quantum version £qs meaning simply to renormalize temperature. The detailed 
calculation, however, shows that this covers only part of the quantum effects and that they also affect the coupling 
between electronic and RC system. In particular, due to the coarse graining (|2f \ the short time expansion (-ft <C 1) 
used in the high temperature limit by Garg et al. is not applicable at lower temperatures where jhf3 ^ 1 . This 
complicates the analysis considerably. 



B. Equations of Motions 



According to the above discussion we restrict ourselves in the sequel to the case x/ = and consider the time 
evolution of P<t s ,cj' (?7 1) on the coarse grained time scale (|2U[1 starting with preparations A(r^) = X(xi = 0, r^). 

Then, the arj-integration in (|I6[) can be performed exactly and one has for the diagonal part of the RC-reduced 
density P(q, t; [x], [r]]) — p(xf — 0, q, t; [x], [r]]) on the coarse grained time scale 



P(q,t;[x],M) = -^~-r / dnX(ri) exp 



1 ~ 

-E ma (g,t,r-i; [x], M) 



(24) 



with the normalization 



m = \ (q^+sar ) (25) 



and the effective action 



1 ^ i ± r i r n r i 1 ( S ( l ) 

-E ma (g,t,r i; [x],[?7]) = -77737 ~ T, Z Z sl7W7 [I 



2 

®r ( M , M ,t)n-$ q ( [x] ,M,t) q -$ ( \x\ ,[v],t). (26) 



Here, the first term describes the equilibrium distribution of the RC with (q 2 ) « l/m[3ui 2 + A the variance of the 
overdamped RC and S(t) is the symmetrized (real) part of the position autocorrelation function 

(q(t)q) = l -{q{t)q + qq(t))-l l -{q{t)q-qq{t)) 

= S(t)+iA(t). (27) 

On the coarse grained time scale these correlation functions simplify to read 

S(t) w ( 9 2 ) cl e" w o 2 */7 A{t) w _ J_ e -o 2 */7 . (28 ) 

2m7 

Eventually, the time dependent coefficients $j ,j = r,q,0 are functionals of the electronic paths and are specified in 
Appendix A. 

The crucial point is now to take into account the impact of the sluggish RC motion onto the electronic dynamics. 
For the spin held fixed, the RC relaxes on one of the diabatic surfaces (see also fig. 1) 

muil ( c 



Vr,(q) = ^ L [q-^2V) (29) 
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(j] = — [+] for donor [acceptor]), the relative position of which is characterized by the reorganization energy 

E T = 2 musk q% with qo = — % . (30) 

A dynamical spin causes transitions between these surfaces that are most likely to occur in the Landau-Zencr (LZ) 
region around q = where the diabatic potential surfaces intersect. For the motion within this region the time scale 
of the RC dynamics is <lz which is the typical time a RC trajectory lives in the LZ range before relaxing towards 
the donor/acceptor minima. The electronic dynamics starts from a diagonal state, the donor, and finally reaches the 
other diagonal state, the acceptor, via jumps through the non-diagonal states ("blips"). The bare spin dynamics can 
thus be visualized as jumps along the edges of a square where the corners represent the spin-Hilbert space. Now, 
this square follows the slow RC motion in such a way that blips appear substantially only if the RC lies in the LZ 
range, while the spin essentially rests in one of its diagonal states away from it. For sufficiently large reorganization 
energies f3E r > 1 the probability for the RC to be activated into the LZ region is small [of the order of the activation 
rate fc ac t ~ ujq exp(— PE T )] and one has a separation of time scales between the dwell times of the spin in its diagonal 
and nondiagonal states, respectively, namely t ac t — 1/fcact 3> ibiip- This scenario allows for a perturbative treatment 
where it suffices to assume h/3 <C tbn P -C £lz < 7/ w o- 

To summarize the situation: We assume an overdamped RC coupled to a spin, activated by a rare event into the LZ 
range where it diffuses for times of order £lz while the spin jumps between diagonal (rj = ±1, x = 0) and nondiagonal 
(r] = 0, x = ±1) states on a time scale of the order of ibn P , i- e -, 

1 7 

-, Tl(3 < tbup < t hZ < -W < tact ■ (31) 

7 

These conditions serve as the starting point for a perturbative approach and will be confirmed self-consistcntly 
afterwards. 

The immediate consequence of l|31fl is that for times t 3> ibiip integrals in the functionals &j ( [ry] , [x] , t) containing 
X-spin paths are approximated by 

dsx(s)f(s)~ ds X (s)f(s)= ds X (t-s)f(t-s), (32) 

Jt-t blip Jo 

where f(s) denotes smooth time dependent functions decaying to zero on the time scale £lz- The typical magnitude 
of these integrals is thus of order tb\i P \f(t)\ and therefore much smaller than integrals of the form 

ds V (s)f( S )^ V (t) f dsf(s) (33) 
Jo 

which are of order ^lz max o<s<t|/( s )l- Of course, the above simplifications do not apply to all possible spin paths, 
but rather only to those that are assumed to give the dominant contributions to the path integral I|1U|) . Accordingly, 
spin dependent terms in the action l|2t)|) of the order of tbiipA or smaller are neglected against terms of order 1 or 
larger (for details see Appendix A). 

Along these lines the action lj2f)|l reduces to 



i ~ T- 1 

-E ma (g,t,r i ; [x],Vf) ' 



2<? 2 > 2foa>(l-f£) 



S(t) ^^o_ 
{q z j moj { 



1 2 

Vf 



2c (<7 2 )ci '" 



+lq ^WJo dSXis) (34) 

where we used rjf = i](t). If we take the limit t — > with the exact S(t) [and not its coarse grained form only 
the first two terms survive: Together with N(t) the second one gives rise to a S function which with 77(0) = 77/ = —1 
restricts to = q+ (co/mcjg) so that one regains from the first term the equilibrium distribution of the RC in the 
donor state. Deviations from this equilibrium are described by the preparation function A(r,). For finite times the 
first two terms describe the dynamics of an overdamped harmonic oscillator in presence of a constant external force 
rjf co/mwg. The last term with the x-electronic path can be interpreted as an effective RC-dependent energy bias of 
the bare electronic system and can thus be incorporated into an effective electronic Hamiltonian 



(35) 
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with k = (q 2 ) / (q 2 ) c \. Accordingly, in the reduced density matrix (|10fl we are now able to re-express the path integrals 
over the electronic paths in terms of matrix elements and obtain for its diagonal part with respect to the RC 



N(t) 



77 #E 



\e* HEh ' ettt \a 



(36) 



where £r,c(<7> t, r^; ijf) comprises the first two (x independent) terms in (|34|l . By replacing A(r^)| — )(— | — > 
Ylcr- a' ^a-i.a'.( r i) \ a i) ( a i\ this result applies also to more general initial electronic states where then A CT()CT <(rj) in- 
cludes also the electronic preparation. 

The final step is now straightforward. We take the time derivative of the matrix elements P a{ , a ' (q, t) on the left and 
the right hand side of l|3tj|) . express them on the right hand side as derivatives with respect to q, and neglect terms 
acting on time scales of order tbn P ■ (ibiip/^Lz) and shorter. This brings us to the central result of this paper, namely, 
the generalization of the ZE to include also the low temperature quantum regime, coined the generalized Zusman 
equations (GZE) henceforth, 



P—(q,t) = £- P_ 



-(q,t) + Y [P-+(q,t) - P+-(q,t)] 

iA 



P++(q,t) = C+ P++(q,t) - — [P-+(q,t) - P+-(q,t)] 
P-+(q,t) = C° P_ + (q, t ) + l -^-[P__(q,t)-P ++ (q,t)]-i 

P+-(q,t) = £° P + .(q,t)- l -^[P^(q,t)-P ++ (q,t)]+i 
Here, the quantum Smoluchowski operators [see 1|22|) ] read 



e + ^ q] P +A q,t). 



rwy dq 



muj o(q - v qo) + ^ 



(37) 



(38) 



and describe the overdamped quantum dynamics of the RC on the donor, the acceptor, and the averaged potential 
surfaces, respectively. The diffusion constant is given by n/ (3 = muj 2 (q' 2 ). Apart from the replacement of the classical 
C c \ by its corresponding quantum operators, the effect of quantum fluctuations also shows up in the Cq/k, dependent 
coupling terms of the non-diagonal matrix elements. The coefficient n contains the equilibrium variance of a damped 
harmonic oscillator 



(q 2 U 



h/3uj 2 



00 dw_ 

-oo ^ 



wy(- 



(cjq — lo 2 ) 2 + lu 2 7 2 (— iu>) 



coth(ujh(3/2) 



(39) 



where 7 denotes the Laplace transform of the classical damping j(t) [see l|f 5|l ] and the integral is determined by the 
real part of the dynamical susceptibility. As discussed in the previous section, for strong damping (|20l) this can be 
written as 



k ~ f + mui P A 



(40) 



with A as in (|23|l . Exploiting the equivalence of reaction coordinate model and spin boson model an alternative 
representation is found to be 



W Jo 



dui Jsb(w) coth(w?i/3/2) 



(41) 



duj J S b(w)/^ 



where Jsb is the spectral density within the spin boson formulation. This can be directly calculated from the spectral 
density I(to) of the RC representation |2lL l25|. For an ohmic spectral density of the form I(uj) = rwyuj exp(— lu/lo c ) 
with a large cut-off frequency uj c 7 one has in the overdamped limit 7/0^0 » 1 a Drudc damping for the spin 
boson system Jsb = ( c o/7) w /( w o/7 2 + uj2 )- Since k > 1 for lower temperatures the effective coupling c^/k between 
electronic and RC system is decreased by quantum fluctuations. Accordingly, the distance between the minima of the 
diabatic surfaces effectively shrinks, thus reflecting the fact that the activation barrier between donor and acceptor 
surface tends to become transparent. Indeed, below in Sec. we will see that k crucially influences the transfer 
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rate at lower temperatures and gives rise to a renormalized reorganization energy E T f k. We note that this is a very 
important point and in contrast to the pop ular approach where quantum effects are taken into account ad hoc by a 
simple renormalization of temperature [2ll |40] • This latter procedure violates the detailed balance condition and is 
apparently not consistent with the derivation of the GZE from the exact dynamics. Of course, in the high temperature 
limit 7^/3 < 1 the GZE Q5) reduce to the ZE. 



C. Range of Validity 

To derive the above GZE we have relied on various assumptions based on time scale separations. While for the 
bare RC dynamics the time scale separation grown out of the overdamped limit is well controlled, for the coupled 
electronic -l-RC motion the conditions (|3 1|) certainly need to be confirmed and also be specified more precisely. 

The time scale £lz is gained from the diffusion constant D of the GZE according to qlz = VDthz where l^z = 
q h~A/E r is the width of the Landau Zener region. With D = n/jm/3 we find i L z = (7 A^o)(A?i/3) 2 / (PkE t ) . To get 
an explicit expression for ibiip is more involved. For A = the GZE for the nondiagonal elements are simply solved 
by calculating its corresponding Greens-function. This has been done for the classical case in j^, |^| so that we 
collect in Appendix B the final results adapted to the quantum case only. The oscillating part (phase) of the Greens 
function is proportional to c$qt which is reminiscent of the time evolution of the nondiagonal projectors, |+)(— | and 
| — )(+|, under the A = O-Hamiltonian coqa z . Accordingly, for the RC fixed (no diffusion) in the vicinity of the LZ 
point q — 0, the corresponding oscillation period i osz tends to infinity for A = and for A / the spin moves 
coherently through its discrete Hilbert space with a typical dwell time of order 1/A in each state. However, when the 
RC diffuses away from the LZ range, t osz decreases. In fact, when the RC approaches the well regions around iqo, 
these oscillations become so fast that on average (averaged over typical times, e.g. of order £lz) the P± T are basically 
washed out and effectively the spin resides in one of its diagonal states. This way, the phase of the nondiagonal Greens 
function encodes the changeover from spin dynamics to spin trapping. Its time scale t osz associated with the diffusive 
RC motion, i.e. |g| w \fiji, is obtained as 

The typical dwell time in a nondiagonal state may consequently be considered as tbiip ~ ^osz- Obviously, for very 
large j/ujq when the RC dynamics almost comes to rest (basically no diffusion), t osz may even exceed 1/A and the 
situation discussed above is regained; the coherent spin dynamics lasts over large periods of times and the typical 
dwell time is given by ibiip ~ 1/A. In this case quantum fluctuations in the RC are negligible and the ZE dynamics 
is valid anyway. Hence, in the sequel we estimate ibiip = ^osz with tbiip <C £lz- 

Let us now come back to our assumptions for the perturbative analysis. The time coarse graining (12211 together 
with the conditions 1|31|1 lead first to the prerequisite for overdamped motion j/u>q 3> 1 and second to 

h 2 6' 2 v A 

« PE r « — , Ahf3 . (43) 

7/^o ^o/7 

By rearranging the terms and taking also into account /3Et > 1 we arrive at the more convenient forms 

E r 1 (TxA . _\ 1/2 A 



Ah/3) 



hA Ah(3 1 \E r ' J uj 2 /j ' 

1 < (3E r « iM-. (44) 
Tip 

Within the overdamped domain 7/c^o ^> 1 these relations define the range of validity for the GZE. As we show in 
Appendix A it is exactly this range where the simplifications of the effective action described above can be applied 
in the sense of a semiclassical perturbation theory. In particular, the first of the above conditions determines a lower 
bound on temperature depending on the ratio E t /%A. Physically, with decreasing T tunneling of the RC becomes 
important that is incorporated in the GZE only in a kind of static approximation (see Sec.0 and the range of validity 
shrinks accordingly, see fig. 3. In comparison the ZE have additionally to obey jhp <C 1 and consequently fail for much 
higher temperatures. The first condition is also readily translated into the spin boson model: for a Debye spectral 
density (Drude model) the full bath cut-off frequency f2 c is identified as fi c w Wq/7 anc ^ ^ ne spin/full-bath coupling 
a as a = E T /2. Then, we observe that the electronic dynamics of the GZE extends from the adiabatic A/f2 c 3> 1 to 
the nonadiabatic limit A/Q, c <C 1 also at lower temperatures (see fig. 3). The second relation is independent of A and 
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gives an upper bound for the reorganization energy. This is understood by recalling that a larger E r leads to steeper 
potential surfaces at the LZ point causing larger momentum transfers. In order for the Smoluchowski reduction to be 
valid, the damping must increase to compensate for that. The same restriction holds in the classical regime and may 
thus explain recent discussions about negative probabilities etc. in the ZE for higher reorganization energies |'23| . In 
particular, it reveals that in a strict sense the GZE/ZE cannot be used for ET in polar solvents where reorganization 
energies are relatively large jU (for a detailed discussion of this point see 

Qualitatively, the electronic dynamics described by the GZE is found to be either coherent or incoherent. The first 
case applies always for sufficiently large A/(u>q/j) ^> 1 when the RC motion tends to be trapped by friction so that 
effectively a static field, distributed according to its equilibrium distribution, is coupled to the spin. As known from 
the classical limit (ZE) this evokes damped coherent oscillations in the electronic occupation probabilities 43] . For 
smaller ratio A/ (0*0/7) and sufficiently f3E r > 1 one reaches the domain of incoherent decay characterized by transfer 
rates. To derive these from the GZE will be the subject of the last section. 



IV. EQUATIONS OF MOTION IN PHASE SPACE 



As we have seen above, for higher reorganization energies the time coarse graining of the RC motion runs into 
conflict with the electronic dynamics. The path integral formulation allows to extract a set of equations in full phase 
space by keeping also nondiagonal elements of the reduced density matrix with respect to the RC (xf ^ 0). Since 
the RC is just a damped harmonic oscillator (coupled to dynamical electronic paths) everything goes through exactly. 
As already mentioned above, for continuous systems this generalization of the quantum Smoluchowski equation to 
a quantum Fokker-Planck equation has already been done in 28] . In particular, one can thus explicitly take into 
account the equilibration process of the momentum. 

Here, starting again with l|16|) and working along the lines described above, we find for the x-independent part of 
the action 



-t, RC {x f ,rf,ri;r)f,t) 



x f r f 



2% 



T 2 

2 X f 



2(q- 



; >( 



1 SW 1 



S(t) 



Co 

2 7 lf 



(45) 



In the low temperature range (p 2 ) ~ (mfvy /~k)\xi{u> c /^i) so that, as expected, nondiagonal elements in the RC are 
strongly suppressed and \xf\ takes only small values, roughly of order 1/^/7. The x-dependent parts of the action are 
approximated as above with the result that to this order of perturbation theory (neglecting terms of order ibiipA or 
smaller) no additional XfX dependent terms have to be retained. 

In a next step we go over to a phase space formulation by taking the Wigner transform of p af ,a' (^Z? r /j 'Oi i- e -> 



'<T f ,<r' f (P, 



1 

2nh 



dx 



■■KP x f 



(46) 



The corresponding equations of motion for the matrix elements of the Wigner density matrix arc then found as 

1A 



(P,q,t) = £ QFP 



*-QFP 



W— 
W„+{p, q,t) 



2 

1A 



W-+(p,q,t)-W + -(p,q,t)] 

(p,q,t) - W ++ (p, q,t)] 
2c 



[W. 



, +-^q) W_ + (p,q,t) 
Tlk 



(47) 



where the equation for W++ follows from the first one by replacing + with — and vice versa, while that for Wl| 

follows from the second one by interchanging + and — and complex conjugation of its last term. 

Compared to the GZE, in the above equations the quantum Smoluchowski operators are replaced by quantum 
Fokker-Planck operators 



^QFP 



d_ 

dp 



[muj 2 (q- r]q ) + jp] 



p d 
to dp 



lip 2 



dp 2 



9 



rn 



d 2 



dpdq 



(48) 



They differ from classical Fokker-Planck operators [38( mainly by the last term that describes coupled pq diffusion 
and has also been found for weak damping |44l |45| . In the high temperature range 771/3 <C 1 its diffusion constant 
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tends to zero due to the equipartition theorem and one recovers from l|47|) . I|48|) the phase space equations derived 
by Garg et al. .21]. For low temperatures jhf3 ^> 1 its impact becomes substantial and is dominated by the growing 
(p 2 ) variance. Numerical results with the quantum Fokker-Planck equation for a simple harmonic oscillator [22| show 
that for strong damping mean values like (p(t)) and correlation functions like (p(t)p) are in excellent agreement with 
exact results and that small deviations only occur in the short time range. Of course, if the momentum is eliminated 
adiabatically from H48J) we regain the GZE. 

With the exact equilibrium variances inserted into (|48|l the corresponding Fokker-Planck operator is even more 
general [45^ : It gives for the set of initial conditions J7J and ohmic friction (supplemented by a high frequency cut-off) 
the exact quantum dynamics of a damped harmonic oscillator in the long time limit independent of temperature and 
damping strength. While in principle the range of validity of (|47|l is defined by the conditions 144|) , one might thus 
be tempted to think that with the "exact" £qfp one could also extend this range. However, a deeper analysis reveals 
that by weakening for lower temperatures one of the conditions in l)44|) additional x an d V dependent terms in the 
effective action [sec (|26|) ] must be taken into account. Equations of motion with time independent coefficients are then 
out of reach. This suggests that the relations l|44(l define necessary and sufficient criteria for the existence of (simple) 
dynamical equations for the coupled electronic/RC motion in the strong damping domain at lower temperatures. 
Work to study the quantum phase space dynamics in detail is in progress and will be presented elsewhere. 



V. QUANTUM TRANSFER RATES 



Electron transfer rates have been calculated analytically and numerically within the RC and the spin boson rep- 
resentation as well. In essence, two limiting cases are distinguishable, in particular since they allow for perturbative 
analytical treatments (cf. fig. 3). Roughly speaking, in the nonadiabatic regime the RC crosses the LZ region fast 
compared to the dynamics of the spin (A/£! c <l,!] c cut-off frequency of the full bath) so that golden rule techniques 
in the electronic couplin g A can be invoked. This way, transfer rates for classical and quantum mechanical baths 
have been obtained [Tsl l46t 147^ . In the opposite case of a very slow RC compared to the bare electronic motion 
(A/£l c ^> 1) the sluggish RC dynamics controls the transfer and established methods like e.g. Kramers flux over 
population approach or transitions state theory on the lower adiabatic surface are applicable P, |24| . A crucial point is 
the changeover between purely adiabatic to nonadiabatic dominated transfer. As first shown by Zusman 20] and then 
later by Garg et al. [2l| the ZE allow to derive an explicit expression for an interpolation formula that reduces to the 
adiabatic (nonadiabatic) rate constant in the respective limits. A detailed survey over the available analytical results 
and their performance compared to real-time Quantum Monte Carlo data was given recently in (2(|: For A/fl c <C 1 
the nonadiabatic rate expression works quite well over the whole temperature range. The same is true for A/f2 c ;>l 
and higher temperatures Ah{3<l where the Zusman rate is valid. The only domain where no explicit rate formula is 
known so far is A/f2 c >l at lower temperatures Ah{3>l. It is exactly this region where the GZE apply. 

The procedure to extract the rate from the ZE has been outlined in [22t |24| . It also works for the GZE so that 
details are omitted here. By solving formally the two equations for the nondiagonal elements (using the proper Greens 
function, see Appendix B) and plugging this result into the two equations for the diagonal ones, we arrive at two 
coupled integro-differential equations. These can be reduced to two coupled partial differential equations by exploiting 
the time scale separation ibiip £lz to end up with 

P—(q,t) = -K(q)[P__(q,t) - P ++ (q,t)}+ C-P—(q,t) 

P++ (q, t) = K(q) [P++ (q, t) - P__ (q, t)} + £+ P++ (q, t) . (49) 
Here, the effective position dependent electronic coupling turns out to be 



A2 poo poo 

K(q) = —Re / dq' / dtG (q', 

1 J-oo JO 



t\q) (50) 



where Go(q',t\q) is the Greens function for the diffusion in the average potential that determines the nondiagonal 
elements. Next, the equations are Laplace transformed (with respect to time) and projection operator techniques 
provide a mapping onto electronic populations. Finally, the transfer rate follows in the long time limit, i.e. for Laplace 
parameter tending to zero, in the form 

fc na 

k+ = r-r r-r (51) 

1 + k^/kf + k rm /k d K ' 

for the forward rate where k™ and k± are so-called nonadiabatic and diffusive rate constants for forward and backward 
transfer, respectively. These are related by the detailed balance condition 

kl/k°:=e f}he (52) 
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with a=na, d. The backward rate fc_ follows from k + by interchanging + «-> — . Specifically, one obtains 

(q±qo) 2 ' 



kT = 



1 Z" 00 
Mr) J-oo 



2(<Z 2 



(53) 



which represents the nonadiabatic rate as an average of the effective electronic coupling over the equilibrium distri- 
bution of the donor and acceptor surfaces, respectively. The diffusion in these harmonic wells is described by the rate 
constants 



i = J dt ^G T {q ,t\q ) V2tt(<7 2 ) exp 



(q±qo) 2 



2(q 2 



- 1 



(54) 



where G T (qo, t\q ) denote the Greens functions for the RC motion in the donor and acceptor (see Appendix B). Both 
expressions can be simplified further for a sufficiently large reorganization energy [3E T > 1. It then turns out that the 
nonadiabatic rates coincide with golden rule rates, i.e., 



kT = 



hA 2 



I TT(3 



exp 



where 



E 



4E t /k 



(55) 



(56) 



are the effective activation energies for the forward/backward rates. The diffusive rates are identical to the overdamped 
Kramers rates for escape on the lower adiabatic surface 




exp 



Accordingly, for the important case of a symmetric transfer (e 
rate k = k + + fc_ , takes the handsome form 



(57) 



0) the ultimate rate for equilibration, i.e. the total 



k = 



A 2 



h 2 TT(3 



1 + g U 4 E t /k 



exp (-/3E T I %k) 



with an adiabaticity parameter 



A 



fJ 



hA 



(58) 



(59) 



For g ^> 1 one recovers the adiabatic, for g <C 1 the nonadiabatic rate constant [cf. (|44|l and fig. 3]. The above rate 
expression looks like the classical Marcus/ Zusman result with a renormalized reorganization energy E T — > E T j ' n. Note 
that this simple renormalization only appears if quantum fluctuations in the co-coupling terms of the equations for 
P± T [see (J2ZJ] are properly taken into account. Since n > 1 quantum fluctuations always reduce the effective energy 
barrier which is to be surmounted. Of course, in the high temperature limit (n — > 1) we regain the known result, for 
lower temperature, however, significant deviations are observed, (i) The ratio n grows with decreasing temperature, 
at very low temperatures roughly linearly with [3. As a consequence, keeping all other parameters fixed, g becomes 
larger with lower T meaning that one approaches the range where the transfer is dominated by adiabatic processes 
already for smaller values of A. This behavior of g is also consistent with the larger range of validity of the GZE at 
lower T compared to the classical case. Moreover, it agrees with recent numerical indications ^^|. (ii) At lower T the 
exponent in the rate expression (5E T j ' n tends to become temperature independent in contrast to the classical result. 
In particular, for an ohmic secondary bath I(u>) = mjui exp(—uj/uj c ) with large cut-off frequency we gain from 1391) 
for large friction 7/0^0 S> 1 



1 



7T7 



[V> (1 + ?j/?7/27r) - V (l + 



(60) 



where -0 denotes the logarithmic derivative of the T- function. The corresponding quantum enhancement k//3 is 
depicted in fig. 4 for various values of the damping strength. Since k enters the exponent in l|58|l . even relatively small 
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deviations from the classical behavior k//3 =1/(3 substantially influence the rate. In fact, if formally the limit of very 
low temperatures jh(3 7$> l,ojQhf3 ^> j/o->o is taken [which in a strict sense is out of the range of validity (|44|l ] . k//3 
saturates and the exponent in l|58|l becomes 

4k 8wol n (7/^o) htoo 

which is identical to the overlap of two harmonic ground state wave packets with overdamped variance 
(2h / rwy) hi^ / co ) / tt localized around ^fqo in the donor and acceptor wells, respectively. Hence, nuclear tunneling 
is included in the above rate formula, at least in a nondynamical way. In comparison with precise Monte Carlo results 
(from 26] ) a remarkable agreement is seen over a broad temperature range from low up to high temperatures (fig 5). 
The effect of nuclear tunneling is clearly observable at low T where the classical rates expression is far off. Only at 
very low T and for A/(a>o/7)~l, i.e. close to the boundary of the range of validity (cf. fig. 3), does the quantum 
Zusman rate overestimate the true rate (fig. 6) but is still much better than the classical rate. In this parameter 
range the dynamics of the nuclear tunneling process is strongly affected by dissipation which is known to reduce the 
quantum rate substantially. To complete this discussion, we note that according to the numerical results in pfij for 
very low temperatures Afr/3 1 and large A/(u)q/j) 3> 1 (cf. fig. 3) transient relaxation dynamics persists up to very 
long times and the equilibration cannot be described by single-exponential decay. Hence, combining the golden rule 
expression [TH E?) with the above quantum Zusman rate covers almost the entire parameter range where a relaxation 
rate can be found at all. 



VI. CONCLUSIONS 



In this work we have studied electron transfer between a donor and acceptor site in strongly condensed phase and at 
lower temperatures where quantum effects are essential. Based on an exact path integral representation of the reduced 
density matrix of the electron/reaction coordinate compound we have employed recent results on the overdamped 
limit in the quantum regime to derive equations of motions for its diagonal part in the RC, named GZE. Utilizing 
the Wigner transform this was extended to include also its nondiagonal part and to obtain a description in phase 
space. Accordingly we have given a rigorous derivation of time evolution equations for an overdamped collective 
vibronic degree of freedom in presence of electronic transitions over a wide parameter range. As another important 
result transfer rates have been gained from the GZE that interpolate between adiabatic and nonadiabatic dominated 
relaxation also in the low temperature domain. 

The fundamental feature that makes the exact non-Markovian quantum dynamics effectively to behave Markovian, 
is a time scale separation in the overdamped RC dynamics influencing the electronic motion as well. The main 
findings are the following: (i) The GZE reduce to the ZE for high temperatures, (ii) At lower temperatures quantum 
fluctuations lead to a renormalized temperature and a renormalized RC-spin coupling that combines to a renormalized 
reorganiztion energy and thus conserves the detailed balance condition, (iii) The range of validity of the GZE can 
be specified in detail and comprises adiabatic and nonadiabatic domains also for Ah/3 > 1. In particular, the 
Smoluchowski reduction requires reorganization energies bounded from above by the damping strength, (iv) The 
phase space Fokker Planck equations include explicitly the momentum equilibration of the RC, but it seems to us 
that they do not allow for a substantial extension of the range of validity beyond that of the GZE. (v) A generalization 
of the Marcus/Zusman expression for the electron transfer rate contains a renormalized reorganization energy such 
that in the low temperature range the rate tends to become temperature independent. It thus describes (at least 
partially) nuclear tunneling in agreement with precise numerical date. The golden rule formula (see e.g. |15ll47| ') and 
the quantum Zusman result give now a sufficiently precise description of the incoherent transfer throughout most of 
the relevant range of parameters. 

What we have not discussed here in detail is the dynamics of the GZE and the quantum phase space equations 
based on numerical solutions. This is certainly an interesting subject, particularly for the latter one, and is presently 
under study. As discussed above, the quantum Smoluchowski limit has also been analyzed for anharmonic potential 
fields so that GZE for these cases may also be found |4^. In any case, the appealing feature of equations of motions is 
the straightforward procedure to solve them. Typically they are not plagued seriously by numerical instabilities and 
can therefore be used to describe the ET dynamics also for very long times. In contrast, while e.g. real time quantum 
Monte Carlo methods provide a numerically exact scheme, quantum noise in form of the so-called dynamical sign 
problem sets a severe limit on the simulation time. Further, the GZE give in some aspects a deeper understanding of 
the underlying physics since they allow for approximate analytical results. 
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Appendix A: Effective Action 



The exact action coefficients $j in l|26|l are read off from the results of [29j. On the coarse grained time scale 
%P, 1/7 <C 7/^0 (Smoluchowski limit) they are given by the following approximate expressions 



4\ 



2g s(t) 
(q 2 ) K{3 
■ 4g 7/^ 
1 (q 2 )d (hp) 2 Jo 

2q s(t) '■' 



du7](u)A(2t -u)+i- 



2 7 /u 2 



[q 2 )d m 2 Jo 
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A(t - u) 
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q 2 o s(t) 
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durj(u)A(t — u) 
du I dsx(u)i](s 
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h/3 A(u) 
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A(u) 
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S(s)S(u) 
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2k [s(t) - 6(s - u)] S(s - u) + s(t)S(2t 



u) 



(62) 



with the abbreviation s(t) 



S 2 (t) 



and A(t) and S(t) as in 



If also the time scale separation h/3 <C 



ibiip ^ thz <C tact is applied, these coefficients reduce further. In a first step all leading order terms in tbiip are kept 
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(63) 



(64) 



with A(Q) = -h/2m-f from (J2HJ- To extract in a second step the relevant terms, we estimate the order of magnitudes 
of the various contributions. There, $ r and Q q carry an overall factor pE T /qo, $0 a factor PE X . Specifically, one 
obtains for <I> r : The 77/ term is of order PE r /q , the x term of order £ PE r /q with 



ip 



A 



k h/mu>o 



(65) 



Next The 77/ term is of order PE T /qo; the first x term is of order tbn P / PE r /qo, the second one of order 
£ PE T /qo. Finally, the <E>o contribution: In the first line terms are of order /3-E r £ or smaller, the second line is of order 
PE T , and the third one of order PE r £tb\i P /hp. It is now assumed that pE T £tbn p /hp < 1 so that all terms of this 
and smaller order (e.g. PE T £) are negligible in the sense of a semiclassical approximation. Accordingly, the expression 
(I34|l for the effective action follows. In fact, within the range of validity of the GZE specified in l|44|) this assumption 
is indeed justified, thus proving the self-consistency of the approach. 
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Appendix B: Greens functions 



The Greens function for the averaged potential is obtained from solving 



d_ 

at 



TlK 



G (q,t-t'\q')=0 



(66) 



with the initial condition lim t ^o Go(q,0\q') = S(q — q') and the boundary conditions G(±oo, t — t'\q') = 0. Details of 
the calculation have been reported in 22, 24]. Here, we specify the final result 



G (q,t\q') = 



1 
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(67) 



where 
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LOn V J KiTl 



(68) 



To estimate ibiip (see Sec. IIII C)l we consider times t <C tLZ — 7/^0 an< ^ retain the imaginary part of the exponent in 
leading order given by the imaginary part to c(q',t). 

The Greens functions for the donor/acceptor surfaces are gained from [<9 t — C^] G^ = and read 



G T (q,t-t'\q') 



2 v /7ra(f) 



exp 



q±q - (q±qo)c ° ( * *' 



4a(t) 



(69) 
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FIGURE CAPTIONS 
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FIG. 1: Diabatic potential surfaces for the reaction coordinate. An equilibrium wave packet is excited by a short laser pulse 
from a dark state (thin line) to the donor state (thick line, minimum at ~qo)- From this nonequilibrium preparation a wave 
packet evolves on the coupled donor and acceptor (thick line, minimum at +qo) surfaces where the transfer occurs in the 
Landau- Zener region near q — 0. 

FIG. 2: Smoluchowski range for a harmonic RC according to l|21|l . The dashed line separates the classical domain -yh/3 <C 1 
from the quantum domain yhp S> 1. 

FIG. 3: Range of validity of the Quantum Zusman equations according to the first conditions in 1441 (shaded) for fixed 
Er/hA = 2. The dashed line separates the classical (small Ahfi) from the quantum range (larger Ahf3). 

FIG. 4: Quantum enhancement vs. inverse temperature for various damping strength together with the classical result k = 1. 

FIG. 5: Electron transfer rates vs. temperature for a symmetric system (e = 0) according to the expression l|58|l (solid line), 
the classical result (k = 1, dashed line), and precise Quantum Monte Carlo data (diamonds, from 26]). E T /hA = 10 and 
A/Q c = 2. 

FIG. 6: Electron transfer rates vs. temperature for a symmetric system (e = 0) according to the expression 15811 (solid line), the 
classical result (k = 1, dashed line), and exact Quantum Monte Carlo data (diamonds, from 26]). E T /hA = 10 and A/f2 c = 1. 
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